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04 ' Abstract 

The issue is the expediency of the QCD notions use in the low energy 
region down to the confinement scale, and, in particular, the efficacy of 
■ the QCD invariant coupling a s (Q 2 ) with a minimal analytic modifica- 

tion in this domain. To this goal, we overview a quite recent progress 
in application of the ghost-free Analytic Perturbative Theory approach 
(with no adjustable parameters) for QCD in the region below 1 GeV. 
Among them the Bethe-Salpeter analysis of the meson spectra and spin- 
| dependent (polarization) Bjorken sum rule. 

The impression is that there is a chance for the theoretically consistent 
and numerically correlated description of hadronic events from Zq till 
a few hundred MeV scale by combination of analytic pQCD and some 
explicit non-perturbative contribution in the spirit of duality. 
This is an invitation to the practitioner community for a more coura- 
geous use of ghost-free QCD coupling models for data analysis in the 
| low energy region. 
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1 The pQCD overview 

X . 

QCD effective coupling a s . Common perturbative QCD (pQCD) based upon 
Feynman diagrams starts with power expansion in a s = g 2 s /^ ~ 0.1 — 0.4, the 
strong interaction parameter analogous to the QED fine structure constant. 

In QFT, an important physical notion is an invariant (or effective, or running) 
coupling function a(Q) , first mentioned in the QED context by Dirac (1933). In 
the current practice it was introduced in the basic renormalization group papers of 
the mid-50s[l]. 

The one-loop invariant QCD coupling sums up leading order (LO) logs into a 
geometric progression (with the Bethke[2j convention for the (3k coefficients) 

1 + a 8 {n)Po H%r) fa L V ' 1271 
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At high enough energy (small distance), the QCD interaction diminishes a s (Q) ~ 
l/\nQ — > as Q/A — > oc; r A -> 0. This feature is the famous phenomenon 
of Asymptotic Freedom. 

At the same time, eq.(l) obeys unphysical singularity (Landau pole) ~ 1/(Q 2 — 
A 2 ) in the low-energy physical region at \Q\ — A ~ 400 MeV. Transition to the 
2-loop case does not resolve the issue. 

The asymptotic freedom behavior 1/ In Q remains dominant in the 2-loop or 
Next-to-Leading-Order (NLO) case. Here, an explicit expression for a s obtained by 
iterative approximate solving [TJ |3] of differential RG equation can be written down 
in a compact the "denominator form" (as it was recently motivated in |4]) 
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with values (3 {A ± 1) = 0.663 =f 0.053 ; ft (4 ± 1) 
The QCD scale in the MS scheme tS n t) = 
happens to be close to the confinement scale A^ 4±1 ^ - 
Rk ~ 10~ 13 cm. 
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Figure 1: Effective QCD coupling correlating all the data in the range from a few GeV 
up to a few hundred GeV. The solid curves correspond to the 2-loop, NLO case. Taken 
from the Bethke paper j^. 



According to Bethke[2], the 2-loop pQCD approximation (T5J) turns out to be 
sufficient for numerical correlation of several dozen of various experiments. Indeed, 
Fig.l gives the evidence for the two-loop pQCD triumph: the NLO theoretical curve 
describes quite accurately - within the current experimental and theoretical errors 
- all the data in the energy range from 5 up to a few hundred GeV. 
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However, below 5 GeV the correlation is not so persuasive. Moreover, in this 
region the data on Fig.l (as well as in the corresponding PDG [5] plot) are rather 
scanty. The reason is not a shortage of experiments but rather troubles of their 
theoretical analysis. Among the latter - the issue of unphysical singularities, like 
the "Landau pole" . 

As it is well known, the widely-used expressions for effective QCD coupling (like 
eqs.flT]),©; see also eq.(7) in Ref. [2]) and eq.(9.5) in Ref.[5|) suffer from spurious 
singularities in the LE physical region at \Q\ ~ A( 3 ) ~ 400 MeV. This trouble is one 
of the main embarrassments for the data analysis by pQCD theory below a few GeV. 

Unphysical pQCD singularity vs. lattice data. At the same time, numerous lattice 
simulation results [6] [8] testify to the regularity of a s (Q) behavior in the region 
below 1 GeV. Indeed, as it was summarized in papers [9j [10], all the lattice data 
indicate smooth growth of a s till specific scale Q = ~ 400 — 500 MeV (that is 
close to A® ) with typical values a s (Q*) ~ 0.5 — 0.8 < 1 . 

This means that common iterative solutions of RG eqs., like (1), (2) not only can 
but should be modified in low-energy region to get rid of singularities and correlate 
with lattice data. 

Modifications of "Common pQCD" in the LE domain. Several attempts to elude 
the pQCD singularities have been undertaken since the 80s. Among them are the 
straightforward freezing [TT] and a few other, more sophisticated, like glueball mech- 
anism [T2] and exponential modification [13J. All of them introduce some model 
parameters. 

Meanwhile, in the mid-90s, an elegant way (free of additional parameters) to 
resolve this issue was proposed by Solovtso^ and collaborators [H] - [16] on the 
basis of the causality principle implemented in the form of the Kallen - Lehmann 
analyticity for the QCD coupling a s (Q 2 ) . Then, on the ground of Q 2 -analyticity, a 
consistent scheme known as Analytic Perturbation Theory (=APT) has been elab- 
orated [17] - [19] during the last decade. 

Below, we give resume of the APT essence (Sect. 2) and its application to data 
(Sect. 3) in the above-mentioned troublesome region. These results rise hopes that 
the Bethke's issue of two-loop a s adequacy can be proliferated to one more order 
of magnitude - down to a few hundred MeV with the help of analytically modified 
QCD coupling and some additional nonperturbative means in the spirit of duality. 

2 Analytic Perturbation Theory 

Here, we start with a sketch of APT. For details see the review papers [2U] - [24] . 
1 Prof. Igor Solovtsov deceased on July 28, 2007. 
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2.1 APT - General 



As it is well known, the 1st step of improving a renormalized PT result is supplied by 
the RG Method [1] which allows one to reveal the correct structure of the singularity 
of a partial solution; in the QFT case - the correct UV and IR asymptotics. Its 
essence is a technique of restoring the so-called renormalizatioi^ invariance. In 
QFT, the RG-improved results obey a drawback, the unphysical singularity. 

In the latter case, the 2nd step, a further improving of RG- invariant PT solution 
should be used. Its main idea, imposing of the analyticity imperative that in turn 
stems out of the causality condition, was first formulated in the QED context |27j. 
A more elaborate QCD counterpart, the APT algorithm, is based on the following 
principles : 

• Causality, that results in the analyticity of the effective coupling in the complex 
Q 2 plane a la Kallen-Lehmann representation^] 

a s (Q 2 ) - a E (Q 2 ) = - I™ do P ^ , . 

This property provides the absence of spurious singularities. 

• Correspondence with perturbative RG-improved input by proper defining 
pip) = Ima s (-(j) . 

• Representation invariance, i.e., compatibility with linear integral transforma- 
tions, like a transition from the Euclidean, transfer momentum, picture to the 
Minkowskian, c.o.m. energy, one: 



POO 
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or the Fourier transition from cxe{Q) to its Distance image q;tj (r) ) that yields 
non-power functional expansions for observables - see, below eqs.(jlJ,(I3J. 



2.2 The APT Algorithm 

Euclidean functions. Euclidean ghost-free expansion functions [28] are defined 

A n {Q 2 ) = r da, p n (a) = - lm[a s (-a - ie)f (3) 

Jo cr + Q 7T 



2 Or, more exactly, by the reparameterization invariance |25j of a partial solution. Recently, this 
RG technique has been devised for a class of boundary value problems of classical mathematical 
physics 

3 For some cases it is implemented in a form of the Jost-Lehmann (see Sec. 4 in Ref . |20j ) repre- 
sentation. 
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via powers of a s . They form a nonpower set of functions {Ak(Q 2 )} that serves as 
a basis for modified non-power APT expansion of RG invariant objects in the Q 
picture, like the Adler D-function. The first of these functions can be treated as an 
Euclidean APT coupling a^Q 2 ) = Ai{Q 2 ) . In the one-loop case 



4\Q 2 ) = i 

Po 



A 5 



ln(Q 2 /A 2 ) A 2 -Q 2 



it differs from the usual one a s (Q 2 ) by the term ~ 1/(A 2 — Q 2 ) that subtracts the 
singularity. 

Here, higher expansion functions are related by the elegant recurrent relation 



1 dA [ n\Q 
nflo d InQ 2 



Minkowskian expansion functions are connected [HJ [T71 |2B] with the Eu- 
clidean ones by contour integral and the reverse 'Adler transformation" 
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A k {-z) ■ A k (Q 2 ) = Q 5 
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The Minkowskian APT coupling om(s) = in the 1-loop case 



ct\ r (s) = — — arccos = 
M V 7 vr/3 Vi 2 + vr 2 



L>0 



1 7T 

— — arctan — , L = ln(s/A 2 ) 
ttPo L 



is quantitatively close to the Euclidean APT one; see Fig. 2. 
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Figure 2: Comparison of singular a s coupling with Euclidean o<e &nd 
Minkowskian otu in a few GeV region. Taken from paper \18y . 
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Non-power APT - Loop and RS Stability. In APT, instead of universal 
power-in-a s (a s (Q 2 ) or a s (s) ) series 

d P t(Q 2 /s) = d 1 a s (Q 2 /s) + d 2 a 2 + 0(a 3 s ) , 

one should use for each representation its own particular non-power expansion 

dUQ 2 ) = di a E (Q 2 ) + d 2 A 2 (Q 2 ) + d 3 A 3 (Q 2 ) + . . . , (4) 

tv(s)= d iaM (s) + d 2 2l 2 (s) + d 3 Vl 3 {s) + ... (5) 

that provides better loop convergence and practical RS independence of observables. 
The 3rd terms in (3), (4) contribute to observables less than 5 % |21j. Again the 
2-loop (NLO) level is practically sufficient. 

Fractional APT. In computation of higher-order corrections to inclusive and 
exclusive processes one deals with non-integer fractional powers of QCD coupling. 
In such a case, special fractional generalization has been devised[29] and successively 
applied to pion form factor [30] and to the Higgs boson decay into a bb pair [3T] . 

2.3 APT functions at LE region 

Comparison of APT Euclidean ole and Minkowskian ctu couplings reveals that 
below 2-3 GeV scale they, being close to each other, differ seriously from the common 
singular a s - see Fig.2. Qualitatively, the same is true for higher expansion functions. 

The APT RenormScheme- and loop- stability. In Fig. 3, we give Euclidean 
APT coupling in the one-, two- and three-loop (NNLO) approximations taken in the 
MS scheme. 

1.5 
1.0 

0.5 
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0.0 0.2 0.4 0.6 0.8 1.0 

Figure 3: Euclidean APT coupling cue in the 1-, 2- and 3-loop cases for the 
MS scheme. Taken from paper JJgf . 

A beautiful feature of these curves is their relative loop stability. The two-loop 
curve below 1 GeV differs from the three-loop one by less than 3 per cent. Hence, 
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the APT two-loop (NLO) curve is accurate enough for practical use at three-flavor 
region. This correlates with the Bethke's thesis for higher energies. 

In a real QCD case, one has to take into account the proper conjunction of 
regions with different values of effective flavour number n f . This, a bit subtle issue 
was elaborated in jT9]. From the practical point of view, one should use common 
matching conditions for recalculation of values at the quark threshold crossing. 
Resulting Euclidean functions Ak turn out to be smooth in the threshold vicinity, 
while Minkowskian ones 21^ remaining continuous have jumps in derivatives. 

Recall here that all this is valid for simple APT functions without additional 
parameters. This version is known as a minimal APT. Below, we shortly mention 
its massive generalization which contains an additional fitting parameter. 

The "massive" APT modification. A quite natural ansatz has been added to 
the minimal APT formalism in [52] ■ There, the lower limit in the Kallen-Lehmann 
integral Eq. (|3]) was changed from zero to m 2 > . This parameter reminiscent of 
pion mass squared can be used for the data fitting (as in Fig.4). 

3 Low energy APT application 

The APT approach during the decade of it existence has been applied to a number of 
low energy (above 1-2 GeV) hadronic observables. One has to mention here sum rules 
[331 El], e + e~ inclusive hadron annihilation [35J, r [36] and T decays[37j, above- 
mentioned for mf actors [30, 31], and some others. For details one could address to 
review papers [21] [22]. Below, we shortly overview quite fresh APT applications to 
processes in a rather low-energy region < 1 GeV. 

3.1 APT and bosonic spectrum 

APT + Bethe-Salpeter formalism. Here, is a summary of recent analysis [38] 
of the meson spectrum by combination of the Bethe-Salpeter equation for the (q, q) 
system with the APT approach. 

By use of the 3-dimensional reduction, the BS eq. takes the form of an eigenvalue 
equation for a squared bound state mass 

M 2 = M 2 + U OG E + f/conf , 

with Mq = y^m 2 + k 2 + \J m\ + k 2 - kinematic term, Uconi ~~ confining potential, 
Uoge - one-gluon exchange potential ~ QCD coupling 

(k|f/ GE|k'> = « S (Q 2 ) M OG e(Q = k - k' , k) . 
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For a given bound state a , one has (for details see Refs.[40j) 

m 2 a = (0 a |M Q 2 |0a) + (0a|^OG E |0a) + (0a | ^Conf 1 0a) • 

Last two relations allow one to extract a s (Q 2 .) values for a low enough momentum 
transfer region 100 MeV < Q a < 1 GeV . 




Figure 4: Comparison of ce s from meson spectrum (points with error bars) and 3-loop 

(3) 

ctE at _ 3 = (417 db 42) MeV (3 solid curves). Singular 3-loop a s coupling (dot- 
dashed) is excluded by data. Dashed lines correspond to the massive APT version 
J32f with m ~ 40 MeV . Taken from paper 1381. 

Results of a s extraction from the bosonic spectrum are given in Fig. 4. 
One sees that meson spectrum data roughly follow a bunch of three curves 
for A^ =3 = (417 ± 42) MeV corresponding to the 2006 world average a s (M|) = 
0.1189±0.0010 . There is also a slight hint at the tendency for BS-extracted devalues 
at Q < 200 MeV to diminish in the IR limit. The dashed curves on Fig.4 just relate 
to this possibility. However, in our opinion, such a scenario is supported only by 
data from the D and F orbital excitations of the (q, q) system. They have big error 
bars and some of them are subject to a doubtful interpretation. If we exclude higher 
states and limit ourselves to the S and P ones, the resulting picture will change. 

APT vs S and P data. In Fig. 5, we show the picture without higher orbital D and 
F excitations. This limited set of data with small error bars quite nicely follows the 
minimal APT coupling curve with the world average 

A^ =3 = 417 MeV value. 

3.2 Bjorken sum rule 

Fresh 2006 Jefferson Lab data on the Bjorken Sum Rule for the moment of the spin- 
dependent structure function r p _ ra i at 0.1 < Q 2 < 3 GeV 2 were analysed recently 
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in the NLO approximation [3T] . Higher twist (HT) values extracted within the APT 
provide evidence for better convergence of HT series as compared to the standard 
pQCD. As a final result, a reasonable quantitative description of the data down to 
350 MeV was achieved. 

Together with the meson spectra evidence FigJHJ this produces an impression 
that minimal APT allows one to enlarge the domain of analytic perturbative QCD 
(supported by transparent non-perturbative elements) description of hadronic events 
down to a few hundred MeV. 

4 APT in QCD: Conclusion 

Meson spectrum data analyzed by the Milano BS-technique with the one-gluon 
exchange potential and confinement ansatz result |39j in a possibility to extract the 
QCD coupling a s (Q) values in the LE domain of momentum transfer Q < 1 GeV. 
In a recent research it was shown [38J HO] that the use of ghost-free analytic QCD 
Euclidean coupling in this analysis yields rather an intriguing correlation (shown 
in Figs. 4 and 5) of the "meson spectrum a s values" in the region 250 MeV < Q < 
1 GeV with the world average a s (M|) . 

Along with this, the arena for the APT nonpower expansion results for the 
Bjorken sum rule [31] is also ranging down as far as to the ~ 300 — 400 MeV scale. 

Both the results - 
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• exclude common a s singular behavior and smooth "freezing" below 1 GeV, 



• support minimal APT extension of pQCD, giving hope for a quasi-perturbative 
consistent quantitative picture from 200 GeV to 200-300 MeV. 

Due to this, there appears a chance for the real possibility of consistent theoreti- 
cal analysis of hadronic processes in the low-energy region, the chance that is based 
on two elements: 

- the procedure of getting rid of spurious singularities, by some low-energy modifi- 
cation of pQGD, like the APT one; 

- addition of some appropriate nonperturbative elements in the spirit of parton- 
hadron duality, like confinement ansatz and higher twist contribution. 

We appeal to the QCD practicing community for a more regular use of ghost-free 
QCD coupling models for data analysis in the low-energy region below 1 -2 GeV. 
Just in this region theoretical errors quite often exceed the experimental ones. 
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